Mott-Hubbard transition in infinite dimensions 
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We analyze the unanalytical structure of metal-insulator transition (MIT) in infinite dimensions. 
A multiple- valued structure in Green's function and other thermodynamical quantities with respect 
to the interaction strength U are observed at low temperatures by introducing a transformation into 
the dynamical mean-field equation of Hubbard model,. A complete description of stable, metastable 
and unstable phases is established in the regime U c i{T) < U < U&iT). The Maxwell construction 
is performed to evaluate the MIT line U*(T). We show how the first-order MIT at U*(T) for T > 
evolves into the second-order one at C/ C 2(0) for T = 0. The phase diagram near MIT is presented. 
PACS numbers: 71.30. +h, 71.27.+a, 71.28.+d 



I. INTRODUCTION 

The Mott-Hubbard metal-insulator transition (MIT) 
is one of the classic topics in strongly-correlated electron 



systems 



1,2,3 



The one-band Hubbard model is thought 



to be the minimum one such that the main features of 
the MIT are grasped. Early studies by Hubbard, and 
Brinkman and Rice explain very well the high and low 
energy behaviors of the local single-particle spectrum, 
respectively. 4 In recent years, the transition between a 
Fermi-liquid phase and an insulating phase is extensively 
studied in the framework of dynamical mean-field the- 
ory (DMFT), 5,6 by utilizing many methods such as iter- 
ative perturbation theory (IPT), 7 quantum Monte Carlo 
(QMC) calculation, 8 " 10 exact diagonalization (ED), 11 ' 12 
projective self-consistent technique (PSCT), 13 and nu- 
merical renormalization group method. 14 ' 15 Below a crit- 
ical temperature T c , 8 a paramagnetic (PM) metallic 
phase and a PM insulating phase 16 coexist in the regime 
U cl {T) < U < U c2 {T), 10 ' 15 ' 17 where U is the on-site 
Coulomb repulsion, and U C \{T) and U C 2(T) are bound- 
aries of the coexistence regime at temperature T. When 
U increases, the metallic phase characterized by a finite 
density of spectrum at the Fermi-level disappears discon- 
tinuously at U C 2(T). On the other hand, the insulating 
phase in the large U regime is destroyed abruptly when 
U decreases to U cl {T). 15 Between U cl {T) and U c2 {T) lies 
a first-order MIT line U*(T), at which the free energies 
of the two coexisting phases are equal. 10 ' 17 At the criti- 
cal temperature T c , both U C \{T) and U C 2{T) equal to U c 
and a second-order transition occurs. Above T c , when U 
increases, the system changes from metallic into insulat- 
ing phase through crossover. 17,18 At zero temperature, it 
is found that the metallic phase has a lower energy in 
the coexisting regime. The MIT occurs at U = U e 2(0) in 
such a way that the energy difference between the two co- 
existing solutions disappears quadratically, i.e., Eg — E^ 1 



- (f/ c2 (o)-c/) 2 . 13 

Despite the enormous efforts on this problem, there is 
no complete description for the stable, metastable and 
unstable phases. These phases always appear near the 
first-order phase transition and are important for un- 
derstanding the behaviors of the system under external 
influences. Traditionally, the first-order nature of MIT 
is disclosed from the hysteresis or discontinuity in some 
relevant physical quantities with respect to U. In this 
way, the unanalytical feature of these quantities is not 
displayed, and the Maxwell construction cannot be per- 
formed explicitly to produce the MIT line U* (T) . Hence 
the MIT deserves further study in this direction. In 
this paper, a transformation is introduced to the DMFT 
self-consistent equation to investigate the unanalytical 
structure in the Green's function, local density of states 
(DOS) and double occupation probability in terms of 
U. For each U in the regime U cl (T) < U < U c2 (T), 
three solutions arise: stable, metastable and unstable 
non-Fermi-liquid metallic solutions. The Maxwell con- 
struction is performed to produce the MIT line U*(T). 
The second-order MIT at T = T c is found to persist at 
a metastable level in the low temperature limit. When 
temperature approaches zero, the first-order transition at 
U*(T) evolves into second-order one 19 at U C 2(0). Finally 
an ED phase diagram of MIT is presented. 



II. MODEL AND METHOD 

We start with the single-band Hubbard model at half- 
filling for the PM solution: 



<i,j>cr i 



(i) 



The notations are conventional. A semicircular bare den- 
sity of states is used: D(e) = (2/-kW 2 )VW 2 — e 2 , and 



1 



W = 1.0 is taken as the energy unit. To introduce our 
method, denote a thermodynamical quantity by Q. For a 
specific temperature, the formal dependence of Q on the 
on-site interaction U is written as Q = Jq(U). We con- 
sider the analogy between MIT and conventional liquid- 
gas transition, which was put forward by Castcllani et 
al.. 20 In this analogy, U corresponds to the pressure P 
as the driving force of phase transition, and D = (n^ni) 
corresponds to the inverse density v. This analogy is 
helpful for us to understand the characters of MIT. 8 In 
the mean-field treatment of MIT, unanalyticity of ther- 
modynamical quantities versus U should appear in a sim- 
ilar way as that appears in the isotherm v = v(P) of van 
der Waals equation. In this way, we realize that for an 
appropriate quantity Q, the function Jq{U) should be 
continuous, but may have "Z"- or "S"' -shaped structure. 
Now we focus on the Green's function at an imaginary 
time t. As an example, we choose Q = G({3/2). A discon- 
tinuous jump in G{(3/2) was observed when we directly 
swept along U axis. In order to avoid such discontinu- 
ity, instead of directly calculating G{(3/2) = fG(p/2)(U) 
for each U, we try to find a solution of G((3/2) for the 
transformed self-consistent equation, 

G(/?/2) = f GW 2)(U X(A G(/3/2))), (2) 

where the parameters A and A are assigned such that 
G(/3/2) is single-valued with respect to U even if the 
function /g(/?/2)(^0 nas a u Z n - or "S"' -shaped struc- 
ture. Existence of such an unanalytical structure is 
strongly indicated by the discontinuity of G((3/2). Af- 
ter the new equation is solved, the original functional 
dependence G(/3/2) = /g(/3/2)(U') is recovered by plot- 
ting the self-consistent solution G((3/2) versus the argu- 
ment V = U — X(A — G{(3/2)). In practice, Eq. (2) 
is combined with the DMFT self-consistent equations by 
iterative calculation. For fixed U, A and A, in each iter- 
ation, we first calculate V = U — \{A — G(/3/2)), where 
G(/3/2) is taken from previous iteration. Take the chem- 
ical potential /i = U' /2 to realize the constraint of half- 
filling. Then U' is used to replace U in the conventional 
algorithm 6 to proceed with the calculation in this itera- 
tion (or from initialization). At the end, besides the new 
set of parameters for the effective single-impurity Ander- 
son model, a new G((5/2) is also produced. Both of them 
are used in next iteration. U' varies with iteration until 
the Green's function converges. Then both Eq. (2) and 
the DMFT equations will be fulfilled. No additional com- 
putational effort is needed in implementing our scheme. 
After the convergence is reached, other quantities, such 
as G(t) for a general r, the double occupancy (n-\n\), 
and the local DOS at the Fermi surface p(0), can be cal- 
culated from the converged Green's function G(iu> n ) and 
V. The whole functional dependence Q = ') can 
be obtained by sweeping U in Eq. (2). This modified 
self-consistency scheme was used in the study of phase 
separation in double-exchange systems. 21 This method 
is also effective to reveal the unanalyticity in MIT, and 



the final results are essentially independent of the selec- 
tion of A and A if only | A| is large enough to "stretch" the 
curve. In this paper, the finite-temperature ED technique 
is used to implement the above scheme. 22 This technique, 
originally proposed by Caffarel and Krauth, 11 uses x 2 fit 
to obtain new parameters of the effective Anderson im- 
purity model in the iteration process. 

III. RESULTS AND DISCUSSION 

The minus imaginary-time Green's function at r = (3/2 
is plotted versus U in Fig. 1(a) for T = 0.01, which is 
lower than the critical temperature T c . It is shown that 
—G(/3/2) = —fG(p/2){U) is indeed continuous with a 
"Z" -shaped structure. Three solutions of G((3/2) coexist 
in an extent of U. At T = 0.01, the results evaluated 
by ED technique of five and six sites are consistent very 
well. We find that our results are independent of suf- 
ficiently large A, the value of A, or the initial seeding 
of iteration. We plot G(r) versus G{(3/2) as shown in 
the inset of Fig. 1(a). The analytical behaviors of these 
curves appear clearly that for each r € [0,/3], G(t) is 
unanalytical in terms of U . Moreover, the coexistence 
boundaries of G(r) have the same values irrespective of 
r. The Fourier transformation G(iu) n ) as well as all other 
quantities will have the same boundary of coexistence, if 
only they depend on the Green's functions analytically. 
So, the boundaries of coexistence U c i(T) and U C 2(T) are 
well defined in MIT. 
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FIG. 1. (a) ^/-dependence of minus imaginary-time 
Green's function — G(t) at r = (3/2, obtained using Ns = 5 
(square)and Ns — 6 (dot), (b) The three coexisting solutions 
of Green's function on Matsubara frequency axis, obtained 
using Ns = 6 at U = 2.47. Both figures are for T = 0.01, 
and the thin lines are for guiding eyes. Inset of (a): — G(t) 
versus — G(/3/2) as U varies (Ns = 6). From top to bottom, 
t = 0/32, /3/16, /J/8, /3/4. 

In Fig. 1(b), three coexisting solutions of the Matsub- 
ara Green's function at T = 0.01 and U = 2.47 are 
presented. The upper and the lower curves are similar 
with those from QMC. 8 ' 10 One is metallic-like (squares) 
and the other is insulating- like (up triangles). A new 
metallic-like solution (dots) is found between them. At 
finite temperatures, this newly discovered solution has 
the highest free energy among the threes and is unstable. 
In this paper, we do not intend to quest for the highest 
numerical precision, but focus on the qualitative features 
of the Mott-Hubbard transition and its temperature evo- 
lution. All our results in the following are obtained using 
Ns = 6 (Ns is the number of sites). 

The minus imaginary-time Green's function —G(t) at 
t = (3/2, the double occupancy D — (n^n{) and the local 
DOS at Fermi surface p(0) are plotted versus U for sev- 
eral temperatures in Fig. 2(a), (b) and (c), respectively. 
In Fig. 2(a), —G((3/2) decreases monotonously when U 
increases at T — 0.04. As the temperature decreases, 
it decays more rapidly in the intermediate U regime. 
At T = 0.025, a singular point arises at U = 2.34 
where the slope diverges. Below this temperature, the 
curve is still continuous, but has a 11 Z" -shaped struc- 
ture. With decreasing temperature, the curve is com- 
pressed along —G((3/2) axis and the coexisting regime 
extends along the U axis. This leads to a more pro- 
nounced multiple-valued structure. We do not find slow- 
ing down of this tendency as temperature is lowered down 
to T = 0.0025. Similar temperature evolution behav- 
iors are also observed in D — U (Fig. 2(b)) and p(0) — U 
(Fig. 2(c)) curves. Within our numerical precision, the 
three kinds of curves produce the same value of critical 
point (U c , T c ) and the same boundaries of coexistence 
regimes for T < T c . This is consistent with our con- 
clusion that the unanalytical structure is universal for all 
thermodynamical quantities. The critical temperature T c 
and interaction U c are thus estimated to be T c w 0.025, 
U c w 2.34, which agree quite well with that obtained 
from QMC: T c = 0.026 ±0.003, U c = 2.38±0.02. 8 At the 
critical point, all three quantities have divergent slopes 
with respective to U. From Fig.2(a)-(c), the boundaries 
of coexistence regime U C \(T) and U C 2(T) are easily ob- 
tained, in contrast to previous approaches. 10 ' 15 ' 

The double occupancy D is of special thermodynam- 
ical significance. The free energy can be evaluated by 
integrating along the D — U lines in Fig.2(b): 

F(U, T) = F(0, T) + f D(U',T)dU'. (3) 
Jo 



In Fig. 2(b), D — U curves have similar unanalytical be- 
havior with v — P isotherms of van der Waals equation. 
At T ^ T c , the curves of D — U are similar with those 
from IPT 23 and QMC. 8 The Mott critical point (U c , T c ) 
was studied in detail by QMC 8 and Landau theory of 
phase transition. 23 At < T < Tc, three solutions co- 
exist in a regime around U c . We compare free ener- 
gies of the three solutions for a fixed U and find that 
the phase with intermediate D has the highest free en- 
ergy. For the metallic (with largest D) and the insu- 
lating (with smallest D) phases, their free energies cross 
at the point U = U*(T), i.e., F M (U* ,T) = Fj(t7*,T). 
We solved this equation numerically to determine U* (T) . 
As U passes by U* from below, a stable metallic phase 
(F M (U,T) < Fi(U,T)) transits into a stable insulating 
phase (F/(J7,T) < F M {U,T)) (See Fig.3). This transi- 
tion is accompanied with a finite jump AD of the double 
occupancy D, which is determined by the Maxwell con- 
struction line as shown in Fig. 2(b). Hence the MIT at 
finite temperatures is identified as a generic first-order 
phase transition. At T < T c , it is interesting that there 
is a discontinuous jump in the slope of the D — U curve 
at U = U cl (T). At this point, dF/dU = D is continu- 
ous but d 2 F/dU 2 = dD/dU is discontinuous. It means 
that the transition is of second order. This singularity 
is directly evolved from that of the critical point (U c , T c ) 
as temperature decreases. It turns out that the second- 
order MIT at T = T c docs not disappear in the regime 
T < T c , but persists to the absolute zero temperature 
at a mctastable level. This feature of MIT is not unique 
among the first-order phase transitions. In the double- 
exchange model for manganites, similar feature was ob- 
served in the isotherm of charge density n versus chemi- 
cal potential \x (See Fig. 2 in Rcf. 21). There, when the 
phase separation (also a typical first-order phase tran- 
sition) between paramagnetic (PM) and ferromagnetic 
(FM) phases appears at T < T c , the second-order PM- 
FM transition persists at the metastable level down to 
zero temperature. 

Now we discuss the phenomenon of critical slowing 
down. In a previous QMC study of MIT, 10 the criti- 
cal slowing down in the convergence of iteration was ob- 
served near the boundaries of coexistence regime, and 
it was used as the indicator to determine these bound- 
aries. Such critical slowing down arises as one tries to 
directly calculate the "Z" -shaped curve for the first-order 
phase transition. In our calculations we did not observe 
the critical slowing down of the same kind since what 
we practically calculated is a "stretched" single-valued 
function /q(J7 — X(A — G(/3/2)), and the coexistence of 
solutions is removed by the transformation. In contrast, 
we observed a weak slowing down of convergence at the 
singular point at U = U C \(T) where the slope of D — U 
curve is discontinuous. This slowing down is associated 
with the metastable second-order phase transition iden- 
tified above. Hence it is different from that discussed in 
Rcf. 10. 
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FIG. 2. (a) [/-dependence of -G(/3/2) at T = 0.04, 0.025, 
0.015, and 0.005. (From top to bottom on the left side) (b) 
and (c) double occupancy D ( in (b) ) and DOS at Fermi 
surface p(0) ( in (c) ) versus U at T = 0.04 (squares), 0.025 
(down triangles), 0.015 (dots), and 0.005 (up triangles). Thin 
lines are for guiding eyes. The thick vertical line in (b) shows 
Maxwell construction for T — 0.005. In small U regime, the 
minor deviation of p(0) at T = 0.005 from 2/nW originates 
from the finite-size effect. 

At zero temperature, Rozenberg 12 and Moeller 13 found 
a second-order MIT at U c — U C 2{0). Here, we discuss how 
the first-order MIT at finite temperature evolves into the 
second-order one at zero temperature. Fig. 2(b) shows 
that when temperature decreases, the linear behavior 
of D — U curve becomes dominant and the intermedi- 
ate branch of D goes closer to the insulating solution. 
As a result, the first-order transition point U*{T) moves 
towards U C 2(T) according to the Maxwell construction. 
This tendency continues when the lowest temperature 
T = 0.0025 in this paper is reached. It is anticipated 
that at T = 0, the intermediate branch of D merges with 
the lowest branch and they become degenerate. The lin- 
ear metallic D — U line meets the insulating branch at 
U c iifi) = U*(0), where the jump of D disappears. Hence 
the transition is of second order at T = 0, and it is con- 
sistent with the D — U curves shown in Fig. 34 of Ref. 
6. The quadratic relation E\ - Ef - (C/ c2 (0) - Uf 



(Ref. 12) is obtained by assuming a strict linear D — U 
curve at T = 0. Here it is worth pointing out that the 
term second- order transition in this context means that 
at C/*(0), dE/dU is continuous and dE 2 /d 2 U is discon- 
tinuous, where E = mm(Eg,Eg 4 ). This transition is 
special in that metastable state still exists near the tran- 
sition point at zero temperature. 

In Fig.2(c), the local DOS at the Fermi surface p(0) 
is calculated by extrapolating —ImG(iuj n )/iT to the limit 
w n — > + . In the small U regime, p(0) does not change 
with interaction and remains its non-interacting value 
2/(nW), as required by Luttinger theorem for momen- 
tum independent self-energy. 5 At temperatures above T c 
(e.g., the curve for T = 0.04 in Fig. 2(c)), the Fermi- 
liquid phase in small U regime evolves continuously 
into an insulating phase in large U regime through a 
crossover. At temperatures below T c , the behavior of 
p(0) in the small and large U regime are similar with 
previous results 6,15 . In the coexisting regime, the metal- 
lic phase at U C 2 is smoothly connected to the insulating 
phase at U c \ through a non-Fermi-liquid phase. Accord- 
ing to the integral along the D — U lines, this non- Fermi- 
liquid phase has the highest free energy. In light of the 
present result, the discontinuity and hysteresis of p(0) 
in Ref. 15 originates from numerical instabilities at the 
boundaries of the coexisting regime. Physically, such in- 
stabilities reflect the unstable character of the metastable 
phases near phase boundaries. The metastable second- 
order transition at U C \(T) is now identified as the tran- 
sition between metallic (p(0) > 0) and insulating (p(0) 
« 0) phases. When temperature decreases, the interme- 
diate solution of p(0) moves downwards, but there is no 
obvious tendency that it will merge with the insulating 
solution of p(0) in the limit T — > 0. This is in contrast to 
the D—U curve as shown in Fig. 2(b). From Fig. 2(c), we 
conclude that in the low temperature limit, as U passes 
by U C 2(0) from below, D varies continuously but p(0) will 
jump to zero from a finite value as was observed in NRG 
study. 14 Other quantities characterizing the MIT such as 
the quasi-particle weight Z (Ref. 11) and local spin-spin 
correlation (M z (t)M z (0)} (Ref. 17) should have similar 
unanalyticity in terms of U. 

Based on the analysis of the Green's function and 
other physical quantities, we can determine the phase 
diagram near the MIT. In Fig. 3, two boundaries of the 
coexistence regime meet at a critical point U c ss 2.34, 
T c w 0.025. Above this point, the metallic phase in 
small U regime evolves into an insulating phase in large 
U regime through a crossover area. Following Ref. 15, 
the U regime where —dp(0)/dU is larger than its half 
maximum is regarded as the crossover regime. Below T c , 
between the two boundary lines three phases coexist, i.e., 
stable, metastable and unstable phases. Positions and 
properties of these phases are labeled in Fig. 3. The 
MIT line U*(T) from Maxwell construction resembles 
that obtained by J. Joo et al.. w Compared with the phase 
diagram of other authors, the U C 2{T) line agrees pretty 
well with NRG results 15 , and the U C \{T) line lies between 
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that from NRG 15 and QMC 10 studies. Although the low- 
est temperature in the present paper is T = 0.0025, our 
phase boundaries extrapolated to zero temperature is in 
good consistency with the results of NRG (Ref.14) and 
PSCT. 13 
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FIG. 3. Phase diagram of Mott-Hubbard transition ob- 
tained from ED method. The squares and crosses mark out 
the coexistence boundaries and the MIT line U*(T), respec- 
tively. Thin lines are eye-guiding lines. Corsses with horizon- 
tal dashed lines above T c show the crossover regime. 

PSCT (9) 

(circle) and NRG (10) (diamonds) results for T = are also 
shown. PM, PI denote paramagnetic metal and paramag- 
netic insulator, respectively. The properties of phases in area 
I, II and I+II are noted in the figure. 



IV. SUMMARY 

In summary, we introduce a transformation into the 
DMFT self-consistent equations to study the unanalyti- 
cal behavior of thermodynamical quantities in the Mott- 
Hubbard metal-insulator transition in infinite dimen- 
sions. Using the ED technique at finite temperatures, the 
"Z" -shaped multiple- valued structure of several quanti- 
ties in terms of U is obtained. An unstable non-Fermi- 
liquid phase as well as the two phases discovered pre- 
viously are found in the regime U C \(T) < U < U C 2(T). 
The MIT line U*(T) is obtained by Maxwell construc- 
tion. The second-order MIT at T = T c persists at the 
metastable level down to zero temperature, and the first- 
order MIT at U*(T) evolves into second-order one at 
C/ C 2(0) when temperature approaches zero. The method 
used in this paper should be useful for the study of other 
first-order phase transitions. 
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